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We consider the minimal two-band model for the Fe-based superconductors with a phenomenological pairing 
interaction which mimics short-range antiferromagnetic (AFM) fluctuations. Two superconducting (SC) gap 
solutions are found to exist with the model: sign-changing s-wave gap (±s-wave) and double d-wave gap states. 
Both solutions hold the approximate relation AJ"'"A^ ~ Af ax N e , a generic feature of two band model with 
a dominant interband pairing interaction. We carried out the calculations of the SC properties of the both 
SC states such as the density of states, temperature dependencies of spin-lattice relaxation rate l/7\, Knight 
shift, and penetration depth, particularly taking into account of the interband coherence factors. The results are 
discussed in comparison with the currently available experimental data. 

PACS numbers: 74.20,74.20-z,74.50 
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I. INTRODUCTION 

Recent discovery of the Fe-based superconducting com- 
pounds provided a great impetus to the research of supercon- 
ductivity (SC). Since the first report on the superconducting 
transition at 7K with the doped LaOFeP by Kamihara et al. 1 , 
various substitutions (mainly, P by As and La by Ce, Gd, Sm, 
Pr) for this mother compound increase the superconducting 
temperature T c over 50K with Sm(OF)FeAs 2 . Intensive inves- 
tigations by many experimentalists and theorists have already 
revealed main metallic and superconducting properties of this 
group of materials. 

From band calculations 3 ' 4 i 5 ' 6 i 7 ' 8 , it is agreed on that the 3d 
electrons of Fe atoms are the main contributors to the conduc- 
tion bands crossing the Fermi surface (FS). Besides the degree 
of degeneracy, the key feature of conduction bands is that it 
consist of hole band(s) around F point and electron band(s) 
around M point [in the notation of the folded Brillouin-zone 
(BZ) scheme^]. 

Regarding the pairing symmetry, there are already many ex- 
periments: (1) Knight shift below T c shows a clear drop in- 
dicating a spin singlet pairing 9 ; (2) tunneling spectroscopy 10 
showed the zero-bias conductance peak (ZBCP) - signature 
of a sign changing gap, but the interpretation of the shape 
of density of states (DOS) iV(co) is diverse; (3) nuclear spin- 
lattice relaxation rate 1/Ti 9 ' 11 unanimously showed no coher- 
ent peak and ~ T 3 dependence below T c , hence strongly sug- 
gesting a d-wave type gap; and (4) specific-heat coefficient 
C(T)/T below T c (Ref. 12 ) - although the measurement is not 
yet reaching low enough temperature - appears T linear indi- 
cating the gap with lines of node. All these experiments ap- 
pear to be consistent with a d-wave type gap. However, recent 
penetration depth measurements with PrFeAsO, Sm(OF)FeAs 
and Nd(OF)FeAs (Ref. 13 ) strongly suggest a fully opened gap 
at low temperatures indicating a s-wave type pairing symme- 
try. 

Regarding the paring glues, the phonon interaction appears 
unlikely mainly because the electron-phonon coupling is es- 



timated to be very weak (A, < 0.2) 14 . On the other hand, this 
series of materials, without doping, commonly has a spin den- 
sity wave (SDW) transition at around ~ 150 K. When the su- 
perconductivity appears with doping, the SDW correlation is 
expected to remain, albeit the long range order disappears. 

Recent neutron-scattering experiments with La(OF)FeAs 
and Ce(OF)FeAs (Ref. 15 ) directly measured the antiferro- 
magnetic (AFM)-type correlation of the Fe d-electron spin 
moment. The overall phase diagram with doping for 
Ce(OF)FeAs reveals a close correlation with an antiferro- 
magnetism and superconductivity, suggesting the important 
role of magnetic fluctuations as a pairing glue. It also 
shows that the generic phase diagram of these compounds 
shares the universal features with the high-r,- cuprates, Pu- 
115 superconductor 16 , and various heavy fermion supercon- 
ductors; namely, the SC occurs in the neighborhood of the 
magnetic long range order when this magnetic order is sup- 
pressed. In particular, the magnetic order is an AFM type. 
This universal phase diagram is very tantalizing because it ap- 
pears to cover a wide class of unconventional SC materials 
with a range of T c from a few mK to 100 K and suggests that 
the AFM fluctuation is a common thread and its characteristic 
energy roughly scales with the SC TJ&. 

For the Fe-based SC materials, several theoretical mod- 
els were already proposed and most of them started with 
the orbital basis of the Fe 3d-electrons including Hubbard 
U interaction(s) and Hund coupling(s) 7 5 - 17 ' 18 ' 19 - 20 . Some of 
these studie s 5,17,18 found that the ±s-wave gap as a domi- 
nant instability. A d-wave gap also often appears as a sec- 
ond instability 1719 . In this paper, we took a phenomenolog- 
ical approach to investigate possible pairing states in the Fe- 
pnitide superconductors. The non-interacting part of Hamil- 
tonian is constructed by choosing a minimal set of topologi- 
cally distinct two bands and the interaction part of Hamilto- 
nian is assumed from the experimental input li— , simulating 
a short range AFM correlation. By solving the coupled BCS 
gap equations, we found the two SC gap solutions: a sign- 
changing s-wave gap and a double d-wave gap. For the both 
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FIG. 1: (Color online) FSs of £/,(£) (red) and £ e (k) (green) bands and 
two gap solutions considered in the paper: (A) ±s-wave gap and (B) 
double d-wave gap. The width of the FS represents the local DOS 
for each band within (£>afm = 20raeV energy. 



SC states, we carried out the calculations of the SC properties 
such as the DOS, temperature dependencies of the spin-lattice 
relaxation rate 1/71, Knight shift, and penetration depth. We 
particularly take into account of the interband coherence fac- 
tors, unique to the two band model, in these calculations. The 
results are discussed in comparison with the currently avail- 
able experimental data. 



H. MODEL 

We propose a minimal phenomenological two band model 
for the Fe-based superconductors. For the non-interacting 
part of Hamiltonian, we observe that several band structure 
calculation s 3 ! 4 ! 5 ! 6 ' 7 ' 8 of the Fe-pnictide compounds reached 
the consensus that the FS of the doped compounds consists of 
two hole pockets and two electron pockets. To keep the essen- 
tial physics, but avoiding unnecessary complexity, we choose 
only two topologically distinct bands: one hole band around 
r point (0,0) and one electron band around M point (±71, ±7C). 

The main phenomenological assumption of our model is 
the interacting part of Hamiltonian. This pairing interaction 
V(q) is chosen to simulate a short range AFM spin fluctu- 
ations peaking at the ordering wave vector Q = (±7C,±7i). 
This assumption is directly motivated by the experimental ob- 
servations of the AFM correlation in the Fe-pnictide com- 
pounds by the neutron-scattering measurements 15,21 . The 
elastic neutron-scattering experiments for La(F x Oi_ A )FeAs 
and Ce(F x Oi_ r )FeAs by Dai and co-workers 15 showed the 
long-range AFM order of the Fe3d-electron spins for the dop- 
ing range of x = — 0.05. When this long-range AFM order 
disappears beyond the doping around x=0.05, SC appears up 
to the doping range x ~ 0.2 (this is only the limit of the mea- 
sured data). This overall phase diagram appears to be generic 
for the Fe-pnictides, and clearly shows that the AFM corre- 
lation is the dominant magnetic correlation in this group of 
materials. 

More importantly, this AFM correlation is expected to con- 
tinue to exist in the doping range where the SC phase occurs 
albeit becoming a short-rang one. This speculation is sup- 
ported by more recent inelastic neutron-scattering experiment 
in the doped (BaK)Fe2As2 compound (7 t =38K) by Christian- 



son et al^i In this experiment, a clear magnetic -resonance 
peak is observed at the expected position of the AFM correla- 
tion, i.e., at (k,k) momentum. With this series of experiments 
it is clear that the AFM correlation is the dominant magnetic 
correlation in the Fe-pnictides among other competing mag- 
netic correlations such as weak ferromagnetism 3 ' 4 , checker- 
board AFM 6 - 8 , and AFM stripe phase 8 ' 22 , which were the- 
oretically proposed. Finally, the coupling matrix element is 
assumed to be a constant for simplicity. The Hamiltonian is 
written as 



H = ^h{k)hl a h ka + Y^^i k ) e lo e kc 

ka ka 

+ E v ( k > k ') h h h -ki h k'i h -k'i 

kk'U 

+ e v ( k > k ') e if-k[ e k'i e -k'i 

kk'U 

+ E v ( k > k ') h I-\ h -ki e k'[ e -k'i 

kk'U 
+ E V{k,k)e\ f "_ k[ h^h_ kl] (1) 

kk'U 

where h ka and e ka are the electron creation operators on the 
hole and the electron bands, respectively. Zh,e (k) are the dis- 
persions of the hole band and electron bands, respectively, de- 
fined as £/,(&) = ^(cosA:,; + cosk y ) + t^coskxCosky + £ /! and 

S-e(k) = tf (cos k x + cos k y ) +f|cos ycos y +e e . In this paper, 
we choose the band parameters as (0.30,0.24,-0.6) for hole 
band and (1.14,0.74,1.70) for electron band with the notation 

The pairing interaction V(k,k ) is phenomenologically de- 
fined below. It is all repulsive in momentum space and it 
represents a short range AFM spin fluctuations as explained 
above. 



V(k,k ) = V M 



\$-X)-Q\* + , 



(2) 



where k and k' are the two- dimensional momenta on the two 
dimensional BZ and the parameter K controls the magnetic 
correlation length as ^afm = 2ko/k (a is the unit-cell dis- 
tance). This interaction mediates the strongest repulsion when 

two momenta k and k' are spanned by the ordering wave vec- 
tor Q. This condition is better fulfilled when the two momenta 
k and k' reside each other on different bands in the model 
band structure (see Fig.l). As a result, the sign-changing s- 
wave gap can form on each band as already suggested by sev- 
eral paper a 5 ' 17 ! 18 ' 19 . However, this opposite sign gap on the 
hole and electron bands is not limited with the ±s-wave state 
([Fig.l (a)]. Another possibility, which conforms to the lattice 
symmetry, is that each band develops a d-wave gap but with 
71 phase shift between two bands [Fig. 1 (b)]. We call this type 
of gap as double d-wave gap. 



We need to mention that our model did not include the 
screened Coulomb interaction (neither did the other theoret- 
ical investigation ;] 5 ' 17 ! 18 ! 19 ), which certainly exists in the Fe- 
pnitide superconductors as well as in all metals in general. 
The screened Coulomb interaction is traditionally treated as 
"Coulomb pseudo potential" // in the conventional phonon- 
driven SC. However, the reliable estimate for its strength is 
practically impossible because even a small difference in //* 
would cause a large change in T c . In the Fe-pnictides, if we 
are to determine T c theoretically, a quantitative estimate of // 
is necessary. We did not include it in our model interaction, 
first, because we do not know how to reliably estimate it in 
these compounds and, second, because the primary purpose of 
the present paper is not the prediction of the precise T c . Still 
we could investigate its generic effects on the different pair- 
ing symmetries such as ± s-wave and double d-wave gaps; 
for example, how large value of // is necessary to kill the ± 
s-wave pairing for a given strength of the AFM interaction. 
We think that this kind of analysis will dilute the focus of 
the present paper and therefore should be a separate investi- 
gation. We briefly remark, however, on the general effects of 
the screened Coulomb interaction. The screened Coulomb in- 
teraction becomes a short-range interaction in real space and 
therefore weakly momentum dependent in momentum space. 
This type of interaction is almost harmless for the d-wave type 
pairing but extremely detrimental for the s-wave type pairing. 
The pairing solutions in this paper should be considered with 
this point in mind. 

Now we solve the Hamiltonian Eq.(l) using the BCS ap- 
proximation and the two band electrons need two SC order 
parameters (OPs) 

k' 

A e (k) = £V(M') < e k , ,e_ t # T > . (4) 

k' 

After decoupling the interaction terms of Eq.(l) using the 
above OPs, the self-consistent mean field conditions lead to 
the following two coupled gap equations. 

A*(Jt) = (5) 

- Z k > [V hh (k,k')A h (k)x h (k')+V he (k,k)Ae(k)x e (k')}, 
A, (t) = 

- L k > [Veh{k,k')A h (k)x h (k')+V ee {k,k')A e (k)x e (k')}. 

(6) 

where Vi,h(k,k ), Vi w (k,k ), etc are the same interaction de- 
fined in Eq.(2) but the subscripts are written to clarify the 
meaning of V hh (k,k') =V(k h ,k' h ), V he (k,k') =V(k h ,k' e ), etc., 
and kh and k e specify the momentum k located on the hole 
and electron bands, respectively. The pair susceptibilities are 
defined as 
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(7) 



where Ei ie (k) — . /q 2 + A^ e (k), and N(0)h, e are the quasipar- 
ticle excitations and the DOS of the hole and electron bands, 
respectively, and g>afm is the cutoff energy of the pairing po- 
tential V(q). 

When we solve the above gap equations Eq(5) and Eq.(6), 
we numerically restricted the momenta kh e ar >d k h . around 
the FSs of the hole and electron bands within OIafm energy 
range. Therefore, the FS shapes and the local DOS A^(0)/ !jg of 
the realistic bands are faithfully taken into account in our gap 
solutions. Also no restriction on the functional forms of the 
gaps Af, e (k) was imposed except the general symmetry de- 
picted in Fig.l, so that the ^-dependence of the gap functions 
A/, e {k) will follow the characteristics of the bands and pairing 
interaction. 



III. GAP SOLUTIONS 

As explained in Sec. II, the main pairing process with the 
AFM spin-fluctuation mediated interaction V(q) is the inter- 
band pair hopping between the hole and the electron bands, 
in which a pair of electrons (k, —k) on the hole band scatters 
to a pair of electrons (k , — k ) on the electron band and vice 
versa. This process is particularly dominant when the size 
of the FS of each band is much smaller than the size of Q 
vector. Considering only this interband pair process [keeping 
only Vhe and V e h terms in Eqs.(5) and (6)], we observe the fact 
that the pair potential A/,(£) for the hole band electrons is pro- 
vided by the pairs of electrons in the electron band and vice 
versa. The physical consequence of it is that the relative sizes 
of the gaps and DOSs on each band are reversed; namely, if 
N h (0) >N e (0), then \A h (k)\ < \A e (k)\ holds in general. This 
relation holds both for the ±s-wave and for the double d-wave 
gap solutions and affects all superconducting properties such 
as tunneling DOS, Knight shift, l/7\, and penetration depth. 

For all numerical calculations in this paper, we choose the 
parameters K = 0.2k (^afm ~ 10a), (Oafm — 20meV, and 
Vm = 10eV (average interaction < V(q) >= 1.115eV). Our 
choice of band parameters produces Nh{0) — 0.74/eV and 
N e (0) = 0.285/eV so that N h {0)/N e (0) « 2.6. We think that 
these numbers represent the Fe-based SC materials but should 
not be taken too seriously; in particular the pairing strength 
Vm = lOeV is chosen freely for demonstration. 



A. ±s-wave gap 

This solution for the Fe-based SC is already proposed by 
several author s 5 ! 17 ! 18 ' 19 Here we demonstrate that this solution 
is indeed realized with a simple phenomenological interac- 
tion, which mimics an AFM spin fluctuations, on the minimal 
two band model representing the Fe-based SC compounds. 

As we described above, the reversed relation between the 
magnitude of the DOSs and the size of gaps holds more rigor- 
ously for the s-wave case and we suggest an approximate re- 
lation A™'W/, w A™ ax N e (see the Appendix for more detailed 
discussions). This relation is a generic feature of the model. 



0.01 - 



> 0.00- 



-0.01 - 



A« 
AW 



6 (in unit of it/2 ) 



■ NJoJ/NJO) 

- • N h („,)/N ol (0) 

- - N (oi)/N (0) 



* I 4 * 

L/i JVJ 



/ 1 
• \ 



/ \ 

i v 



added, the ground state of our model is the ±s-wave SC state. 
This conclusion is already obtained by other authors 5 ' 17 ' 18 ' 19 
with different models and approaches. The separate and total 
DOSs plotted in Fig. 2(b) show the main features of the its- 
wave gap: two peak structure, the large DOS with a small gap 
and the small DOS with a large gap. The overall shape of the 
total DOS is not very much revealing compared to the current 
tunneling DOS measurements 10 . However, it is too early to 
make a decisive conclusion with our calculations without in- 
cluding Andreev scattering. Also the ZBCP, the hallmark of 
a d-wave gap and observed in experiments with the Fe-based 
superconductors 10 , can equally be obtained with the ±s-wave 
gap state. 

We consider nuclear spin-lattice relaxation rate l/7\ for 
the ±s-wave gap. Several groups 911 have reported that l/7*i 
shows no coherence peak and the T 3 power law below T c , 
strongly suggesting an unconventional gap with lines of node 
such as a d-wave gap. s-wave gap is known to have a construc- 
tive coherent factors for \/T\ to induce the coherence peak 
over a temperature range below T c . However, as Mazin et al. 5 
envisaged, the sign-changing gaps between two bands provide 
a destructive coherent factor for the interband scattering which 
will largely cancel the intraband coherent factors. As a result 
the coherent peak of 1 /T\ for the ±s-wave gap SC will be sub- 
stantially reduced. The explicit formula that we used for the 
calculations is the following: 



m (in unit of A max ) 



FIG. 2: (Color online) (a) ±s-wave gap solutions A/,(£) and A e (k). 
(b) Normalized DOS of the hole band N/, (to) (red dotted line), elec- 
tron band N e (co) (blue dotted line), and the total N, or (co) (solid black 
squares). 



Therefore, given a substantial difference of DOS between the 
hole and electron bands (several band calculation ;; 3 ' 4 ' 5 ' 6 ' 7 ' 8 in- 
dicate that this is true for the Fe-based SC materials), at least 
two distinctively different sizes of the SC gaps should be ob- 
served in various experiments 23 . In particular, because the 
band with a larger DOS would dominate the physical proper- 
ties but actually holds a smaller gap, this feature will modify 
various SC properties of the Fe-based SC in unorthodox man- 
ner, such as A/T c value, temperature dependencies of various 
SC properties below T c , and the responses to impurities. 

In Fig. 2, the solution of the ±s-wave gap and the corre- 
sponding DOS are shown. As mentioned, 2V/,(0) = 2.6N e (Q) 
for our bands. Accordingly the size of gaps of the hole 
band and the electron band is reversed as A™' 1 s=s lOmeV and 
A™ ax m 25meV. The strongly momentum dependent pair- 
ing interaction and the realistic bands naturally induce an 
anisotropic modulation of the s-wave gaps with C\ symmetry; 
the modulation is stronger for the larger gap on the smaller 
DOS band (electron band around M point). Compared to the 
case of the double d-wave solution, the average size of the 
±s-wave gap is larger by a factor of ~ 5 with the same pair- 
ing potential. Therefore, unless some other interactions are 
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For the temperature dependence of the gaps Ah ie (k,T), we 
use a phenomenological formula, A/,, e (k,T) = Ai ie (k,T = 
0) tnnh($y/T c /T-l). By choosing the values ofA%f/T c , we 
can partially take into account of the strong coupling super- 
conductivity effect, p is not a sensitive parameter for final 
results; we take p = 1 .74 in this paper. 

Fig. 3(a) shows the contributions to the l/7i relaxation rate 
from each terms of Eq.(8): the hole band, the electron band, 
and the interband terms. It shows that the cancellation of the 
coherence factors is not perfect in general unless the condi- 
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FIG. 3: (Color online) 1/7] (T) of the ±s-wave gap with A™ av /T c 
= 1.5. (a) Separate term contributions of Eq.(8): total (solid black 
square), hole band (open red square), electron band (open blue tri- 
angle), and interband term (solid green square), (b) Total l/7i(T) 
without (solid black square) and with (open green pentagon) damp- 
ing, (c) The log-log plot of (b). The inset is a wide view 



tions N h (0) =N e {0) as well as \A h (k)\ = \A e (k)\ are fulfilled. 
Nevertheless, due to the large cancellation by the interband 
coherence factor, the height of the coherence peak is very 
much reduced [compare the total l/T\ and the hole band only 
l/T\ in Fig. 3(a)]. Small amount of impurities can easily wash 
out this reduced coherence peak as shown in Fig. 3(b); the 
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FIG. 4: (Color online) Normalized Knight shift (uniform spin sus- 
ceptibility) of ±s-wave gap. The total (solid black square), hole band 
(open red circle), and electron band (open blue triangle) contribu- 
tions are shown separately, (a) Kj™* /T c = 1.5 and (b)A™ ax /T c =0.5. 



damping rate T = 0.05A™"* is enough to completely kill the 
coherence peak. The subtle part is to fit the ~ T 3 power law 
below T c , It requires to tune R — A™' v /r c ratio. In Fig. 3, 
R = 1.5 (automatically, it makes A'"" x /T c « 3.75 which is 
quite a large value) is used for the best fit. Fig. 3(c) shows 
that this pseudo-T 3 behavior is not extended to the very low- 
temperature region as in the d-wave case because this T 3 be- 
havior in the ±s-wave gap is not an intrinsic property of the 
lines of nodes. At low temperatures, there appears the ex- 
ponential drop inevitably due to the full gaps, and then it fi- 
nally reaches the impurity-induced T-linear region because 
we added some amount of impurities to kill the coherent peak. 
All these details put rather stringent conditions to confirm the 
±s-wave gap state with experiments. 

Now we consider the Knight shift which is the measure of 
uniform susceptibility in SC phase. Because it is a q — ► 
probe, there is no interband contribution and the total Knight 
shift is just sum of the contributions from each band as fol- 
lows: 




FIG. 5: (Color online) Normalized superfluidity density 
X (0)/A, (T) of ±s-wave gap and its separate contributions 
from the hole and electron bands, (a) A^ x /T c = 1.5 and (b)A^ OT /r c 
=0.5. 
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In Fig. 4, the normalized Knight shift (uniform spin sus- 
ceptibility) is plotted and it shows the typical flat behavior 
of a s-wave gap at low temperatures. The contributions from 
the hole and electron bands show separately the feature of 
the larger DOS with small gap and the smaller DOS with a 
larger gap. Fig. 4(A) is the results with A™ ax /T c =1.5, the 
best gap-T c ratio to fit T 3 behavior of \/T\ in Fig. 3. Fig4(B) 
showed the results with A'™ x /T c =0.5 to demonstrate the con- 
vex shape of Knight shift which was reported by Matano et 
ali^ for Pr(FO)FeAs. Impurities does not change much of this 
feature unlike in the case of d-wave gap. 

Finally, we consider the penetration depth. The static 
response function to the electromagnetic fields is the 
following: 24 




The q = limit of this kernel K(q — 0, T) is directly pro- 
portional to the superfluidity density or \I\\(T) in the Lon- 
don limit. For our two band model, total kernel is the sum of 
K/,(q,T) and K e (q,T) with the proper weighting factor pro- 
portional to the DOS N/, e (0) of each band and there is no 
interband screening current, a = (vf/2)qk is the non-local 
parameter and can be rewritten in more convenient form as 



and the penetration depth at zero temperature, respectively. 
Apparently if ($2-) = OCo is small compared to 1, the non-local 
effect becomes negligible. A typical value of OCo for YBCO 
was estimated about 0.01, for example 24 . For the Fe -based 
superconductors, we believe that oco is not much larger than 
the values of the high-r c cuprates. Also for a s-wave case, 
the non-local effect does not change much of the temperature 
dependence of 1/A, 2 (r) except the overall magnitude. There- 
fore, we take (Xo = 0.0 for the calculations of the penetration 
depth in the ±s-wave case. However, this effect can induce 
an important modifications in the d-wave case, which will be 
discussed in Sec.IIIB. 

Fig. 5 shows the normalized superfluidity density 
A. 2 (0)/A. 2 (r) and separate contributions from the hole 
and electron bands. The exponentially flat region appears 
at low temperatures due to the full gap opening, which is 
consistent with recent experiments. 13 Relatively narrower 
region of the flat part (for T < 0.271) compared to the 
ordinary s-wave gap is another feature due to the smaller gap 
with the larger DOS of the ±s-wave gap SC. A subtle part 
here is to fit the high- temperature region (0.3 T c < T < T c ). 
With A'™ x /T c =1.5 (the same value used for the \/T\ fit), 
this part becomes too convex [Fig. 5(A)] in comparison to the 
experiments. 13 A smaller gap-!T 6 ratio can make it concave as 
shown in Fig. 5(B) (with A'™ x /T c =0.5); this concave feature 
was recently observed by Martin et al^ 

In summary, the ±s-wave gap state provides the most con- 
sistent descriptions for the penetration depth experiments. 
However, it explains l/7\ only for a limited temperature 
range, even with a fine tuning of A/ he /T c ratio and impurities. 
Knight shift of any shape can be fit with two band parameter 
(this is also true with the double d-wave gap). The tunneling 
DOS does not provide a decisive conclusion. 



B. Double d-wave gap 

In Fig. 6, the gap solution and the corresponding DOS of 
the double d-wave gap are shown. As mentioned, our model 
bands have A^,(0) = 2.6N e (Q), and consequently gap in the 
2meV is smaller than the one of the elec- 
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vf /A' nax and A-o are the coherence length 



h 

tron band A" ,ax rs 4-meV. The sizes of the maximum gaps are 
~ 5 times smaller than the ±s-wave gap solutions. There- 
fore, the double d-wave gap solution is not the best SC state 
for our phenomenological model with an antiferromagnetic 
pairing interaction. This result is in agreement with other the- 
oretical studies 5, iiJl. This conclusion may change with the 
correlation length of the AFM fluctuations, the sizes of the FS 
of the hole and electron bands, etc. But we numerically found 
that ±s-wave gap solution is favored compared to the double 
d-wave gap solution for most cases. As discussed in Sec. II, 
however, the screened Coulomb interaction may change this 
conclusion. The detailed studies about this issue will be dis- 
cussed in a separate paper. 

To complete the comparisons, we calculated the same SC 
properties of the double d-wave gap state. The separate and 
total DOSs shown in Fig. 6(b) demonstrates the main features 
of the double d-wave gap: the large DOS band with a small 
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FIG. 6: (Color online) (a) Double d-wave gap solutions A/,(£) and 
A e (k). (b) Normalized DOS of the hole band N/,(co) (red dotted line), 
electron band N e (co) (blue dotted line), and the total N /0 , (co) (black 
squares). 



FIG. 7: (Color online) 1/71 (T) of the double d-wave gap. (a) Each 
term contributions of Eq.(8) : total (solid black square), hole band 
(open red square), electron band (open blue triangle), and interband 
term (open green triangle). (b)Log-log plot of total 1/71 (T) for 
A™*/r c = 1.75 (black square) and 1.0 (green circle). 



gap and the small DOS band with a large gap. This result 
shows a similar feature of the tunneling DOS measurement 
by Wang et al. 10 except the ZBCP, which does not show up in 
our simple DOS calculation but should appear when the tun- 
neling conductivity is properly calculated with Andreev scat- 
tering process. 

We consider the nuclear spin-lattice relaxation rate l/7i. 
As in the case of ±s-wave gap, there are three contribu- 
tions for the total 1/71 relaxation rate: hole band, electron 
band, and interband scattering terms. The formula is the 
same as Eq.(8) but the last three terms should be dropped 
because the FS average of A/, e (k) vanishes in this case. We 
use the same form of temperature depe ndent gap function as 
A h ,e(k,T) = A h}e (k,T = 0)t<mh(1.74^T c /T-l). In Fig.7, 
R = A™77>1.75 is used for the best T 3 fit below T c . How- 
ever, in the double d-wave gap, R = 1.5—2.5 provide reason- 
ably good fits, showing a more tolerance than the ±s-wave 
gap state. 

Fig. 7(a) shows the separate contributions from each chan- 
nel together with the total contribution. As in the ±s-wave 
case, \jT\.h provides the largest contribution and \/T\ e pro- 
vides the smallest contribution. There is also the interband 



term 1/T/' 



In contrast to the ±s-wave case, all three 



terms display a similar temperature dependence and no coher- 
ence peaks. Fig. 7(b) shows the same \/T[" tal (black squares) 
in log-log plot. The overall features of 1/7] are the ones 
of the typical d-wave SC state: no coherence peak near T c 
and ~ T 3 below T c , and consistent with the current NMR 
experiments 911 . At very low temperatures, T-linear behav- 
ior starts to appear due to a small damping for the numerical 
calculations (r = 0.005A™"). For comparison, we also show 
Yjjtotai ( p en g reen circles) with A"™ x /T c =1.0, a smaller 
gap-T c ratio; it exhibits a substantial convex part below T c and 
then starts displaying the universal T 3 behavior before enter- 
ing the impurity dominating region. 

In Fig. 8, we show the result of the uniform spin suscep- 
tibility which is measured as Knight shift. Fig. 8(a) shows 
the results with A"™ x /T c = 1.75. The hole band contribu- 
tion is dominant as in 1/71 and the electron band contribu- 
tion show the steeper drop just below T c because of the larger 
gap-r c ratio A'" ax /T c w 3.5. The overall behavior of the to- 
tal % S (T) below T c shows a typical d-wave behavior such as 
T-linear at low temperatures. Fig. 8(b) shows the results with 
A™ /T c = 1.0. A smaller gap-T c ratio makes the Knight shift 
convex as in the ±s-wave case and observed by Matano et al. 9 
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FIG. 8: (Color online) Normalized Knight shift (uniform spin sus- 
ceptibility) of the double d-wave gap: the total (solid black square), 
the hole band (open red circle), and electron band (open blue trian- 
gle) contributions are shown separately, (a) A™ ax /T c = 1.75 and (b) 

A™7r e =i.o 



for Pr(FO)FeAS. This result demonstrates that this convex be- 
havior of Knight shift is irrelevant to the gap symmetry but a 
generic feature of the two-gap (or multigap) SC. But it reveals 
that the gap-7 c ratio A" WX (Q)/T C is much smaller than the stan- 
dard BCS value, where A mflx (0) refers to the gap of the band 
with largest DOS. 

Now we calculate the penetration depth. As we discussed 
in Sec. II, most of experiments, up to now, report a flat tem- 
perature dependence of A,(7) at low-temperature region and 
suggest a fully gapped SC state 13 . A naive double d-wave gap 
state has no chance to explain this flat behavior at low temper- 
atures. Therefore, we consider a non-local effect of the elec- 
tromagnetic response of the double d-wave gap superconduc- 
tor as a possible cause to modify the typical temperature de- 
pendence. In order to include the effect of the non-local elec- 
trodynamics, we use the fully q— dependent kernel Ki ie (q, 7) 
[Eq.(10)] and put it into the integral formula for A,(7) with the 
specular boundary condition, 
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FIG. 9: (Color online) Normalized superfluidity density 
X (0)/X (T) of double d-wave gap and its separate contribu- 
tions from the hole and electron bands with A™ aA /T r = 1.75. (a) 
5o 



oc = ^ = 0.0 and (b) a =0.5. 



where Ki re (q,T) are the normalized kernels as Kf, e (Q,0) = 1 
and q = qXo is a dimensionless momentum. The results with a 
diffusive boundary condition are qualitatively the same; there- 
fore, they will not be discussed. For the non-local parameter 
ao, we think it to be much smaller than 1 for the Fe-based su- 
perconductors, but here we take it as a free parameter and see 
how large value is needed to fit experimental data. 

Fig. 9 and Fig. 10 show the normalized total superfluidity 
density A, 2 (0)/A, 2 (r) and separate contributions from the hole 
and electron bands for the double d-wave gap state. Fig. 9 used 
A™7r t =1.75 and Fig.10 used A™ ax /T c =1.0. In each figure, 
panel (A) is a local limit (a = 0.0) and the panel (B) is a non- 
local limit (a = 0.5). The local cases display the typical d- 
wave behavior at low temperatures, i.e., the linear decrease in 



T from T = 0. The extreme non-local cases (ao 



So _ 
*o 



0.5) 



introduce a substantial round-off (~ T 2 ) region at low temper- 
atures, which is, however, not an exponentially flat behavior 
as the recent experiments claim. Further, even a rough fitting 
requires an unreasonably large non-local parameter ao. 

In summary, the double d-wave gap state can provide con- 
sistent descriptions for tunneling DOS, 1/71, and Knight shift. 
However, there is an intrinsic difficulty to explain the flat be- 
havior of the penetration depth at low temperatures. Also, in 
our model with an AFM mediated pairing interaction only, the 
double d-wave gap solution is energetically less favored than 
the ±s-wave gap solution. 



IV. CONCLUSION 

We demonstrated that a minimal model with a phenomeno- 
logical pairing interaction of the AFM spin fluctuations can 
allow both the ±s-wave gap and the double d-wave gap solu- 
tions with the realistic bands of the Fe-based SC compounds. 
With the same parameters, the ±s-wave gap solution is ener- 
getically more favorable by a factor of ^5 times, so that it has 
a better chance to be realized in the Fe-based SC compounds. 
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FIG. 10: (Color online) Normalized superfluidity density 
X (0)/X (T) of double d-wave gap and their separate contributions 



from the hole and electron bands with A"' ax /T c = 1.0. (a) (Xo = 
0.0; (b) Oto =0.5. 



FIG. A. 1: (Color online) The normalized gap ratio |A C |/|A/,| vs the 
dimensionless coupling constant X = \/Vi ll ,V l ,i 1 Ni 1 N e . The bottom 
baseline is sjNijN e and the top baseline is N[,/N e . 



\max 



In both cases, we found that the approximate relation 
Nh ~ A'" ax N e holds because it is a generic feature of the 
two gap SC when an interband pair scattering is the dominant 
pairing interaction. This relation appears for all SC proper- 
ties in subtle way, which modifies the value of A(0)/7 C and 
other SC properties in unorthodox way. Numerically solving 
the coupled gap equations for the two bands, we found the 
detailed structure of the gap functions Af,, e {k), which showed 
an anisotropy (~ 20 %) of the ±s-wave gaps. We also cal- 
culated the key SC properties, for both gap states, such as 
tunneling DOS, 1/71, Knight shift, and penetration depth and 
discussed them in comparison with experiments. When we 
calculated these quantities, we paid a special attention to the 
interband coherence factor which is a unique feature of multi- 
gap SC. This interband coherence factor particularly produced 
an important modification to the 1/71 relaxation rate of the 
±s-wave gap state. 

The ±s-wave gap state provides the most consistent de- 
scriptions for the penetration depth experiments: the flat low- 
temperature behavior. 13 Besides the low-temperature behav- 
ior, the high-temperature (0.3 T c <T < T c ) behavior - due to 
a large difference of the gap sizes A/, and A e and their corre- 
sponding DOS Nh e - can be either concave or convex. How- 
ever, l/7i experiments can only be fitted for a limited tem- 
perature range, even with a fine tuning of A/, e /7 c ratio and 
impurities. The ±s-wave gap state is not inconsistent with 
the Knight shift and the tunneling DOS data, but overall does 
not provide any decisive merit in comparison with the double 
d-wave gap. 

The double d-wave gap state, although it is energetically 
less favorable in our model unless additional interactions are 
added, provides the best fit to the 1/71 experiments. However, 
it has a difficulty to explain the penetration depth experiments 
for the low temperature flat behavior. It requires an unreason- 
ably large non-local effect to fit the low temperature part; it is 
still not exponentially flat but only ~ 7 2 . If this low temper- 
ature part of X(T) is, indeed, confirmed to be exponentially 



flat, the double d-wave gap state should be ruled out. Tunnel- 
ing DOS and Knight shift can be fit with the double d-wave 
gap state as much as with the ±s-wave gap state. 

In conclusion, quantitative calculations, carried out in this 
paper, with the two most promising SC gap states can serve 
as guidelines for sorting out the possible pairing states of the 
Fe-based SC in comparison with the current and future ex- 
periments. For that, very low temperature measurements and 
systematic studies with the amount of impurities will provide 
decisive information to determine the correct gap symmetry. 

Note added - Recently, we have known that similar stud- 
ies of 1/71 for the ±s-wave state were carried out by two 
groups 25 where only the interband scattering process was ana- 
lyzed and by another recent paper 26 where both the interband 
and intraband processes were considered as in our paper. 
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APPENDIX A: RELATION BETWEEN A/,/A e AND N h /N e 

In the main text, we claimed the approximate relation 
NhAh sw N e A e as a generic feature of the two- band model with 
a dominant interband interaction. This kind of relation will 
have direct and important implications to the experimental ob- 
servations. However it is pointed out by Mazin 27 that this is 
not a rigorous identity in general. In this appendix, we clarify 
the degree of the validity of this relation. 

Here we consider the ±s-wave gap state only. Assuming 
constant gaps, A/,(k) = A/, and A e (k) = —A e , and only the in- 
terband interaction, the coupled gap equations (5) and (6) are 
simplified as 



10 



A/, = V he A e x e {T,A ei (a AFM ) 
A e = V eh A h Xh(TAh,(OAFM) 



(A.l) 
(A.2) 



where %t, and % e are defined with Eq.(7). 

First, when T =T C , the above equations can be written as 



A/, = VheN e A e x const 

A e = VehNhAi, x const 



(A3) 
(A.4) 



C m AFM 



&, 



tanh 



§ 



logl.l4a)AFAf/7c. Because 



where const = J m, — t- 

of Vfe = Vp/,, we immediately obtain the following relation 

from the above equations. 



A e m, 

& III i rji rji 

— = \ / — ; when T = T C . 

A h V Ak 



(A.5) 



Next, when 7* = 0, Eqs. (A.l) and (A.2) are written as 



A/, = V he N e A e \og( 



&AFM + 


y^AFM 


+AI 




A e 




®AFM + 


y^AFM 


+ Aj, 



(A.6) 



A, = V eh N h A h log( y — — -) (A.7) 

A/, 

In general, these equations do not yield a simple algebraic 
relation between A e / 'A/, and Nf,/N e , but we can obtain the 
simple relations for the limiting cases. First, for the extreme 
weak coupling limit, ie., when A/, e <C ®afm, the two logarith- 
mic terms become asymptotically equal as \og(2(i>AFM / A e ) « 
log(2cOAFM/A/,), and we obtain the same relation as the T = T c 
case Eq.(A.5). On the other hand, for strong-coupling limit, 
ie., when A/, >e ^> COafm (which is certainly an unphysical 

limit), log([CdAFM + J®afm + A h,e}/ A h,e) ~ ®AFM/A h , e and 
we obtain the relation, 



Having found the results of the two limiting cases, we can 
guess that the gap ratio A e /Aj, should be in between these two 
limiting ratios. For example, we can attempt an expansion 
with x = log [Nh/Ne] starting from the weak coupling limit 27 , 
and we obtain, in the first order of x, 



IN, log [AW] 

Ne { 4 



(A.9) 



where X = yVhe^^tWe is a dimensionless coupling con- 
stant. For practical use, we numerically solve Eqs. (A.l) and 
Eq.(A.2) and plot the ratio A^/A/, as a function of A, in Fig. A.l. 
The ratio A e /A/, becomes a universal curve when it is normal- 
ized by the distance between J -^- and -^-. The result indeed 
shows that when X ss 1, it is in between two limiting ratios, 
jj 1 and jj 1 as we expected from the above analysis. 

In reality, there exist two complications. First, the intraband 
couplings Vhh and V ee need to be included. A little analysis 
of Eqs. (5) and (6) as well as of numerical results reveals 
that this effect always enhances the gap ratio toward the limit 
jj 1 . Another complication arises from the fact that there are 

more than two bands in real material a 3 i 4 ' 5 i 6,7 i 8 . Applying the 
results of the above analysis, we can suggest the following 
approximate relations. First, we classify the bands of the real 
materials into two groups: the hole bands around F point and 
the electron bands around M point, respectively. Then in the 
strong-coupling limit, 



£A /v M,,«^A e ,,^, (A.10) 

i i 

and in the extreme weak-coupling limit, 

I>mV^7~E A ^V^- (A- 11 ) 



A, 
A/, 



Nh_ 
' N, 



when T = and A/,, e > &afm- (A.8) 



Considering several uncertainties in reality, Eq.(A.lO) can 
serve as a practical rule of the thumb. 
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